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The synergy between the Fermi-\jPC£ and ground-based Cherenkov telescope arrays gives us the 
opportunity for the first time to characterize the high energy emission (fOO MeV - fO TeV) from 
more than 30 blazars. In this study we performed a Fermi-hNT spectral analysis for all TeV-detected 
blazars and combined it with archival TeV spectra. Our results for low synchrotron-peaked BL Lacs 
(LBL) show hints of absorption features in the GeV band that could be interpreted as internal 
opacity at the source. We note that simple or broken power laws cannot describe all the observed 
GeV- TeV spectra and more complex spectral shapes seem required. 



I. INTRODUCTION 

Active Galactic Nuclei (AGN) are extreme extra- 
galactic objects with an observed luminosity outshin- 
ing their host galaxy. Their non-thermal continuum 
emission extending from radio band to X rays or 7 
rays suggest underlying relativistic emission mecha- 
nisms. AGN have been subject to comprehensive 
studies since their discovery, leading to classification 
and unification schemes based on multiwavelength ob- 
servations and polarization measurements In this 
frame, AGN are classified as either radio-quiet or 
radio-loud, and blazars constitute a subclass of the 
latter, with their jet axis oriented close to the ob- 
server's line of sight. This particular orientation com- 
bined with relativistic beaming gives rise to prominent 
observational features in blazars, such as anisotropic 
radiation, rapid variability, high polarization and su- 
perluminal motion. Blazars are divided into two sub- 
classes, flat spectrum radio quasars (FSRQ) and BL 
Lacertae objects (BL Lac). FSRQs are characterized 
by a broad line emission region, depicted with strong 
lines in their optical spectra, that are not present in 
BL Lacs. 

The broadband spectral energy distribution (SED) 
of blazars exhibit a two-component structure, with a 
low-energy component peaking around IR to UV band 
and a high-energy one around X-ray to 7-ray energies. 
The underlying emission mechanism reponsible for 
low-energy component is believed to be synchrotron 
emission from relativistic electrons in the blazar jets. 
On the other hand, the high-energy component could 
either have leptonic 0, [J or hadronic 4, 5] origin. 
Instrumental selection effects brought up two distinct 
classes for BL Lacs as radio-selected (RBL) or X-ray- 



selected (XBL). This instrumental classification was 
later replaced with a more physical one based on ra- 
dio to optical and optical to X-ray spectral indices, in- 
troducing the terminology of high- and low-frequency- 
cutoff BL Lac objects (HBLs and LBLs) An- 
other approach based on the peak frequency of the 
synchrotron component of the SED is also used to de- 
fine the same classification (high- and low-frequency- 
peaked BL Lacs as HBLs and LBLs). Detection of 
intermediate objects (IBLs) between these two obser- 
vationally distinct groups has made it more plausible 
that BL Lac objects constitute a continuum rather 
than a discrete sequence (e.g.,[7]). 

Mkn 421 was the first blazar and extragalactic 
object detected as a very high-energy (VHE; E > 
100 GeV) 7-ray emitter with the Whipple telescope 
in 1992 Q. Since then, different candidate selection 
methods have been applied to radio. X-ray or high- 
energy (HE; E > 100 MeV) 7-ray blazar data in the 
aim of finding new TeV blazars [^-[ll|. To date, 40 
blazars have been detected in the TeV sky[25|, with a 
census consisting of 29 HBLs, 4 IBLs, 4 LBLs and 3 
FSRQs. 



II. SOURCE SAMPLE 

Our blazar sample contains all blazars with a pub- 
lished VHE spectrum before February 2011, with a to- 
tal of 26 sources (see Table|T|. This includes 19 HBLs, 
3 IBLs, 2 LBLs and 2 FSRQs. Seven of these blazars 
were detected with EGRET and 23 of them are in the 
Fermi 2-year catalog. More than half of the sample 
have been detected multiple times in the VHE band. 
These multiple detections extending over several years 
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and obtained mostly with different instruments sug- 
gest that spectral variability in the VHE band is a 
common property for VHE blazars. Even though no 
general pattern has been established for VHE vari- 
ability, several sources have been observed to have a 
flux increase up to a few times their baseline emis- 
sion p^ - [l^ . occasionally accompanied by a change 
in spectral index and minute-scale flux doubling 
times [Ulil. 



Name 


SED type z Fermi state 


RGB J0152-f017 


HBL 0.080 average 


3C 66A 


IBL 0.444? simultaneous 


lES 0229+200 


HBL 0.140 average 


lES 0347-121 


HBL 0.188 average 


PKS 0548-322 


HBL 0.069 average 


RGB J0710+591 


HBL 0.125 simultaneous 


S5 0716+714 


LBL 0.300 high 


lES 0806+524 


HBL 0.138 average 


lES 1011+496 


HBL 0.212 high 


lES 1101-232 


HBL 0.186 average 


Markarian 421 


HBL 0.031 average 


Markarian 180 


HBL 0.046 average 


lES 1218+304 


HBL 0.182 simultaneous 


W Comae 


IBL 0.102 high 


4C +21.35 


FSRQ 0.432 simultaneous 


3C 279 


FSRQ 0.536 high 


PKS 1424+240 


IBL 0.260? high 


H 1426+428 


HBL 0.129 average 


PG 1553+113 


HBL 0.4? high 


Markarian 501 


HBL 0.034 low 


lES 1959+650 


HBL 0.048 low 


PKS 2005-489 


HBL 0.071 average 


PKS 2155-304 


HBL 0.117 simultaneous 


BL Lacertae 


LBL 0.069 high 


lES 2344+514 


HBL 0.044 average 


H 2356-309 


HBL 0.167 average 



TABLE I: VHE blazar sample. No secure redshift is avail- 
able for 3C 66A, PKS 1424+240 and PG 1553+113. The 
Fermi states "low" and "high" are as described in the text. 
In some cases simultaneous data were available, and in the 
remaining cases 27-month time-averaged Fermi data were 
used. 

The first 27 month of Fermi data and archival VHE 
spectra published before February 2011 were used to 
construct combined GeV-TeV SEDs. These include 
six data sets where VHE data overlaps with the Fermi 
era (RGB J0710+591, lES 1218+304, PKS 1424+240, 
PKS 2155-304 and two different measurements for 3C 
66 A). The remainder of the VHE data were taken be- 
fore the Fermi mission. All VHE spectra were cor- 
rected for the extragalactic background light absorp- 
tion using fl5| . 



III. FERMI ANALYSIS 

The fact that most of the GeV and TeV data are 
not contemporaneous introduces caveats for the inter- 
pretation of the combined spectra. Moreover, Fermi 
data represents an average state over a fairly long pe- 
riod, whereas the VHE spectra consist of day-scale 
"snapshots", mostly taken during flares. As a solu- 
tion to this problem, for bright-enough sources, the 
Fermi data were split into "low" and "high" states 
as described below. Thus, non-contemporaneous GeV 
and TeV measurements were matched in a more re- 
alistic way than directly using all the time-averaged 
Fermi data. 

Diffuse class events with energy between 300 MeV- 
100 GeV from the first 27-month of Fermi data (from 
4 August 2008 to 4 November 2010) were used for 
the analysis. For blazars that were observed during 
the Fermi era, only the time periods of a few months 
that cover the corresponding VHE observations were 
selected. 

First, a 27-month light-curve analysis was per- 
formed for each source using an aperture photometry 
technique. Events from a region of 1° radius from the 
target location were selected and counts were plot- 
ted as a function of time, each time bin containing 
49 counts. For sources with high statistics, low- and 
high- flux states were identifled and separated using 
the average count rate as a threshold (see Figure [T] 
and Table U). 
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FIG. 1: Aperture photometry light curve of the HBL 
Mrk501. The solid line represents the average count rate, 
that is used as a reference point to split the data into low- 
and high-flux states. 

Next, a spectral analysis was done for each data set. 
Events from a region of interest of 8° were selected and 
analyzed with Fermi Science Tools v9rl8p6^265, us- 
ing instrument response functions P3_V6_DIFFUSE. 
Sources from the IFGL catalog [l^, bright spots 
with test statistics > 25 and standard galactic and 
isotropic diffuse emission background components [27j 
within the region of interest were included in the 
source model flies. Spectral points were calculated 
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using an unbinned maximum-likelihood analysis tech- 
nique [13 • Figure [5] shows the combined SED for 
the HBL MrkSOl. 
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PL: r = 1.95 + 0.02, NDF = 32.7 / 10 

BPL;r, = 1.83±0.04,r =2.71±0.17,E = 266.1 ± 72.7, NDF = 2.5 / 8 
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FIG. 2: Combined GeV-TeV SED of the HBL Mrk 501 
(black circles). The VHE spectrum belongs to a low flux 
state Therefore, the low-state Fermi spectrum was 

selected for the combined GeV-TeV SED. A broken power 
law yields a better fit than a simple power law. 




FIG. 3: Distribution of peak frequency bands. 



He II line (54.4 eV) 



IV. IC COMPONENT 



H I & He II combined 



The peak frequency of the inverse- Compton (IC) 
component is a salient parameter for describing blazar 
non-thermal continua and studying population trends. 
Systematic studies for measuring the IC peak fre- 
quency mostly suffer from the lack of statistics and 
simultaneous data [lO]. In this work, we focus on 
finding the IC "peak frequency band" rather than 
the "peak frequency" , using a model independent ap- 
proach. For each blazar SED, we identify the en- 
ergy decade in which the largest amount of power is 
emitted. Figure [3] shows the distribution of peak fre- 
quency bands. We observe that the low-synchrotron- 
peaked objects (LSPs) have the maximum of their 
emission mostly below 1 GeV. On the other hand, 
high-synchrotron-peaked objects (HSPs) tend to peak 
in the TeV range. 



V. INTERNAL ABSORPTION 

A power-law function often fails to describe the 
combined SEDs successfully. For most blazars, the 
fits are significantly improved with a broken power 
law (see Figure [5]) or an absorbed power law. Partic- 
ularly, some LBL objects display absorption-like fea- 
tures in the 10-100 GeV range. To investigate this 
feature, we considered the following absorption sce- 
narios from the broad line region (BLR) 21] for the 
blazars S5 0716+714 (LBL), W Comae (IBL) and BL 
Lacertae (LBL): 



Full BLR spectrum 



In all the cases, absorption from the He II com- 
plex seems to be dominant. A double-absorption sce- 
nario does not provide an improved fit over the He 
II single-line absorption. On the other hand, the full 
BLR absorption underproduces the TeV emission (see 
Figure |3]). We also modeled a full BLR absorption 
combined with cascade emission from the 7-7 absorp- 
tion [22,]. No significant contribution from cascading 
effects in the Fermi band is expected (see Figure [5]). 



H|+He||:i;„ = 1.7±2.0 t„^=3.5+1.2 
BLR:t = 9.9±5.9 log ^ = 2.5 
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FIG. 4: Fermi high-state + MAGIC Q spectrum for the 
LBL S5 0617+714. The red curve represents H I + He II 
absorption only and the blue one the full BLR absorption. 
The power-law fit is shown with the dashed line. 



H I line (13.6 eV) 
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months of Fermi data and combined it with archival 
TeV data. We see absorption-like features mostly in 
LBL spectra, these featm^es are best described as BLR 
absorption from the He II complex. In addition, we 
observe that the peak frequency in the GeV-TeV re- 
gion follows the synchrotron peak frequency trend. 
Complete results will be discussed in a forthcoming 
paper. 



FIG. 5: Fermi high-state -I- VERITAS [24| spectrum for 
the LBL W Comae, with absorption of a pure power law 
(blue) and a power law with exponential cut-ofF (red), by 
the full BLR emission, with cascading effects taken into 
account. 



VI. SUMMARY 
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In this paper we give some results from our ongoing 
GeV-TeV blazar population studies. We analyzed 27 
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